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Abstract 

We consider two models of random diffusion in random environment in two 
dimensions. The first example is the self-repelling Brownian polymer, this describes 
a diffusion pushed by the negative gradient of its own occupation time measure 
(local time). The second example is a diffusion in a fixed random environment 
given by the curl of massless Gaussian free field. 

In both cases we show that the process is superdiffusive: the variance grows 
faster than linearly with time. We give lower and upper bounds of the order of 
t log log t, respectively, tlogi. We also present computations for an anisotropic ver- 
sion of the self-repelling Brownian polymer where we give lower and upper bounds 
of i(logt) 1 / 2 , respectively, tlogt. The bounds are given in the sense of Laplace 
transforms, the proofs rely on the resolvent method. 

The true order of the variance for these processes is expected to be i(logi) 1 / 2 for 
the isotropic and i(logt) 2//3 for the non-isotropic case. In the appendix we present 
a non-rigorous derivation of these scaling exponents. 
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1 Introduction 

We consider two models of random motion in random environment in d = 2: 

• self-repelling Brownian polymer process, abbreviated in the sequel as SRBP, re- 
spectively, 

• diffusion in the curl of massless Gaussian free field, abbreviated in the sequel as 
DCGF. 
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In both cases the critical dimension of the model-class is d — 2: for d > 3 the displace- 
ments are diffusive (i.e. the variance grows linearly in time) and for d = 2 multiplicative 
logarithmic corrections are expected. We provide a rigorous proof for the superdiffusivity: 
we give lower and upper bounds of order t log log t, respectively, tlogt on the variance 
of the displacement in the sense of Laplace transforms. The lower bounds are the more 
interesting, the upper bound being almost straightforward. 

The SRBP model is continuous space-time counterpart of the so called true self- 
avoiding random walk (TSAW). The class of models has a long history. It first appeared 
in the theoretical physics literature where the models were formulated, and based on 
scaling and renormalization group arguments dimension dependent scaling behavior was 
conjectured. See [3], [H], |15| . Shortly later and seemingly totally independently of 
the physics papers mentioned, similar questions arose in the probabilistic literature. See 
|13) . |1] where further conjectures and partial results appeared. For a concise historical 
account of the models, conjectures and results we refer the reader to the survey [21] and 
the introduction of the more recent paper [6]. Our present result completes the picture 
in the sense that after the d — 1 and d > 3 behaviour being more-or-less clarified in [20j, 
|23] , |22) . |19j . respectively, in [5J, we now settle the question of superdiffusivity in d = 2. 
However, note that our bounds are still far from being sharp. There is plenty of room 
left for improvements. 

The DCGF model belongs to the class of random walks and diffusions in random 
environment. The drift field, being the curl of a scalar field in 2d, is divergence free and 
hence the environment seen by the random walker is a priori stationary. There are too 
many papers on this topics to list here. An instructive and rich survey, though not the 
most recent one, is [10]. For a more recent survey see chapter 11 of [8]. It turns out 
that the model considered here, where the drift field is the curl of 2d massless Gaussian 
free field (mollified by convolving with a smooth approximate delta-function) is just 
on the boundary between diffusive and superdiffusive asymptotics. The more robustly 
superdiffusive cases were considered in [HJ. Our results complete the picture presented in 
[9] in the sense of proving superdiffusivity in this borderline case. 

The expected order of the diffusivity for both of these models is t(logt) 1 / 2 and it 
is conjectured that this is universal among a wide class of two-dimensional isotropic 
models. (See [11] for another example.) In the appendix (Section H]) we present a non- 
rigorous, nevertheless very instructive scaling argument which explains this conjecture. 
The argument dates back to the late sixties, cf. [I], [2], [5J. In our opinion this sheds 
sharp light on the origins of superdiffusivity in tracer diffusion models. The argument 
shows that the isotropy of the model is important: in particular for non- isotropic two- 
dimensional models the argument gives t(logt) 2 / 3 for the diffusivity. This was rigorously 
proved for the diffusivity of a second class particle for finite range asymmetric exclusion 
models by Yau [24J building on the results of [TJ]. To illustrate the difference we also 
consider a non-isotropic version of SRBP, for this model we give lower and upper bounds 
on the diffusivity of order t(logt) 1 / 2 respectively tlogt. 

The structure of the paper is the following. In the remaining parts of the Introduction 
we fix notation (subsection II. ip . formulate the models in rigorous mathematical terms 
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(subsection II. 2p . describe the picture of the environment as seen by the random walker 
(subsection II. 3p and formulate the main result of the paper (subsection ll.4p . In section [2] 
the relevant Fock space and the relevant operators therein are presented (subsection 12.11) 
and the variational problem is explicitly formulated (subsection 12. 2p . Section [3] contains 
the computational parts of the proof: we give upper bounds on the nontrivial term in the 
variational problem (subsection 13. ip and prove the main lemmas needed for completing 
the proof of the main result (subsection I3.2p . Subsection 13.31 contains the computational 
parts of the proof for the non-isotropic model. Finally, section H] is an appendix which 
contains the mentioned non-rigorous scaling argument. 

1.1 Some general notation 

Throughout this paper we are in two dimensions. V : M 2 — > IR will denote a (fixed) 
approximate delta-function, that is a smooth (C°°), spherically symmetric function with 
sufficiently fast decay at infinity (exponential decay certainly suffices). We also impose 
the condition of positive definiteness of V: 



V(p) = [ e ip - x V(x)dx > 0. 
Jr 2 



Occasionally we will also use the notation U : IR 2 — > M for the unique positive definite 
function which yields 

V = U * U, U(p) = \jv{p). 

A particular choice could be 

e -M 2 /(2<x 2 ) e -M> 2 
V(x) = — - — - — , U(x) 
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but the choice of the concrete function V is of no importance. 

Partial derivatives (with respect to space coordinates) will be denoted 9j, % = 1,2. 
The gradient and curl of a smooth scalar field A : M 2 — > R are the vector fields 



gradA = (d 1 A,d 2 A), curl A = (d 2 A, -d x A). 
The divergence and rotation of a smooth vector field A : IR 2 — > IR 2 are the scalar fields 

div A(x) = d 1 A 1 + d 2 A 2 , rot A(x) = d 2 A x - d x A 2 . 
Given a vectorial expression b = (pi, b 2 ) we will denote 

b =(6 lj 6 2 ) = (6 2 ,-6i). 
1 1— >■ B(t) will always denote a two-dimensional standard Brownian motion. 



3 



1.2 The displacement processes considered 
The self-repelling Brownian polymer process (SRBP) 

Let x i — y F(x) be a smooth gradient (that is rotation free) vector field on R 2 with slow 
increase at infinity. We define the stochastic process t i— >• X(t) 6l 2 as the solution of 
the following SDE: 

dX(t) = ^F{X{t))- J gvadV{X{t) - X(u))du\ dt + V2dB{t). (1) 

Introducing the occupation time measure (also called local time in this paper) 

lit, A) = |{0 < s <t : X(s) G A} | 

where A C R 2 is any measurable domain, we can rewrite the previous SDE as follows: 

dX(t)= (F-gradV *£(t,-))(X(t))dt + y/2dB(t). (2) 

The form (|2J) of the driving mechanism, shows explicitly the phenomenological meaning 
of the law of the process: it is pushed by the negative gradient of its own local time 
towards domains less visited in the past. 

The non-isotropic SRBP 

Let now F\ : IR 2 — > R be a smooth scalar field on R 2 with slow increase at infinity. The 
process is again defined as a solution of an SDE similar to (Q]), but instead of grad we 
use the operator (d\, 0): 

dXx{t) = (^Fi(X(t)) - J dxV^Xit) - X(u))du^j dt+V^dB^t), dX 2 {t) = V2dB 2 (t). 

Or, expressed in terms of gradient of local time: 

dX x {t) = (F x - d x V * l(t, -)){X{t))dt + y/2dB x (t), dX 2 (t) = V2dB 2 (t). 

Now, only the first coordinate of the displacement is pushed by the corresponding negative 
gradient component of mollified local time. The second coordinate moves as memoryless 
Brownian motion. 

Diffusion in the curl of Gaussian free field (DCGF) 

Now let x t— > F(x) be a smooth curl (that is divergence free) vector field on R 2 with slow 
increase at infinity. We define the stochastic process t i— > X(t) G R 2 as the solution of 
the SDE 

dX(t) = F(X(t))dt + V2dB(t). (3) 
This is a diffusion in the drift field F. 
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1.3 The environment seen from the moving point 



SRBP 

The environment profile appearing on the right-hand side of (j2J), as seen in a moving 
coordinate frame tied to the current position of the displacement process is t i— > r](t, •): 

rj(t, x) = (F- grad V * £(t, •)) (X(t) + x) 

with initial value 

?7(0, x) = F(x). 

t t— > rj(t, ■) is a Markov process with continuous sample paths in the Frechet space 
tt= {cu E C°°{R 2 M 2 ) : rot a; = 0, || u \\ k ^ r < 00} 



where || u \\ k mr are the seminorms 



sup {l + \x\y 1/r \d^ md u k (x)\ (4) 



x& 



defined for k — 1, 2, multiindices m = (m 1; . . . , m<f), rrtj > 0, and r > 1. 
Let 7r be the Gaussian measure on Q defined by the covariances 



u k (x)d7r(uj) = 0, K k i(x-y)= / a; fe (a;)a;i(2/)d7r(a;) = F * g kl (x - y) 
n Jn 



where 



2 1 X k X[ 



x 
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0*j(aO = -<9 fcZ log|x| 
The Fourier transform of the correlations is 

km = r%v{ P ). 

\p\ 

It is clear that u distributed according to 7r the gradient of the massless Gaussian free 
field smeared out by convolution with U. 

It has been proved in [B] that the Gaussian probability measure it is stationary and 
ergodic for the Markov process rj t . That is: if the initial vector field F is sampled from 
this distribution, then the vector field profile seen from the position of the moving particle 
will have the same distribution at any later time. (Although [6] deals with the d > 3 
case the same proof applies for d = 2. See also [19] for the Id case.) The process will 
be considered in this stationary regime. That is, the initial profile rj = F is distributed 
according to the stationary measure n. 

Non-isotropic SRBP 
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The situation is similar to the previous case, but our environmental profile will be the 
scalar field f : R 2 ->■ R 

C(t,x) = (F - d 1 V * £(t,-))(X(t) + x). 

This will be now a Markov process with continuous sample paths on 

n = {coe C°°(R 2 R) : || u \\ mtr < 00} 

and the seminorms || uj || are defined very similarly to (EE1) . 

Computations very similar to those in [19] and [6] show that the Gaussian measure 
n defined by the expectations and covariances 



ou(x)dir(uj) = 0, K(x — y)= / uj(x)oo(y)dir(ui) = V(x — y) 
a, Jn 

is stationary and ergo die for this Markov process. 

In order to keep unified and consistent notation, in the anisotropic case we will still 
use vectorial notation for the scalar field u, as follows: 

uii{x) = cu(x), u)2(x) = 0. 

In this case the probability measure it will be denoted as Gaussian distribution of a two 
component vector field with covariances 

K M {x -y) = V(x - y)6 ktl 6 ltl . 



DCGF 

The environment seen by the moving point is now t t— > rj(t, •) defined as: 

rj{t,x) = F(X{t)+x). 

Now, 1 1— >■ T)(t,-) is again a Markov process with continuous sample paths in the Frechet 
space 

n = {cue C°°(R 2 R 2 ) : divu = 0, || u < 00} 

where the seminorms || uj |L m _ are formally defined as in (j3J). 

Let now tt be the Gaussian probability measure on Q defined by the expectations and 
covariances 

u k (x)dir(uj) = 0, K k i(x - y) := / u k (x)uji(y)d7r(uj) = V * g ki (x - y), 
n Jn 



where now 



9ki{x) = -9fc,log|x| 



x 



3 • 
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The Fourier transform of the correlations is now 

\v\ 

The random vector field Uki{x) is now the curl of the massless Gaussian free field smeared 
out by convolution with U . 

It is a well known fact that the the process rj(t) of a diffusion fl3]) in any translation 
invariant and ergodic divergence free vector field F is (time) stationary and ergodic (see 
e.g. [TO], or chapter 11 of [8]). In particular, this holds if the drift field F is sampled 
from the distribution tt. 

1.4 Superdiffusive bounds 

For both processes considered we denote 

POO 

E(t) = E (| X(t) | 2 ) , E(\)= E(s)e- Xs ds. 

Jo 

The main result of the present paper is the superdiffusive lower bound stated in the 
following theorem, which is valid for both processes. 

Theorem 1. Let X(t) be either one of the (isotropic) processes SRBP or DCGF started 
with initial environment profile sampled from the respective stationary distributions. There 
exist constants < C\, C% < oo such that for < A < 1 the following bounds hold 

d\~ 2 log | log A | < E(X) < C 2 \~ 2 | log A | . (5) 

Remarks: (1) Modulo Tauberian inversion, these bounds mean in real time 

C 3 tloglogt < E (\X(t) | 2 ) < Citlogt, 

with < C3, C4 < oo and for t sufficiently large. 

(2) Based on the Alder- Wainwright argument sketched in the Appendix (Section H]) the 
expected true order is E (| X(t) | 2 ) x t(logt) 1 / 2 . 

The proof of Theorem [U follows the main lines of |12j . However on the computational 
level there are some notable differences. The bounds exploited in [12] which rely inter alia 
on a clever application of Schwarz's inequality don't yield superdiffusive lower bounds 
in our case. This is due to the isotropy of our model as opposed to the anisotropy of 
the asymmetric simple exclusion models. From the Alder- Wainwright scaling argument 
it follows (see the Appendix) that for anisotropic models E (| X(t) | 2 ) x t(logt) 2//3 , while 
for isotropic models E (| X(t) | 2 ) x t(logt) 1 ^ 2 is naturally expected. Thus, the isotropic 
models are less superdiffusive than the anisotropic ones. This phenomenon manifests also 
on computational level: proof of superdiffusive lower bound for isotropic models exhibits 
some subtle a priori extra difficulty. 

For the non-isotropic model we have the following bounds: 
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Theorem 2. If X(t) is the non-isotropic SRBP then we get the bounds 

Ci\~ 2 1 log A | 1/2 < E(X) < C 2 A~ 2 |logA|. (6) 
Remark: Modulo Tauberian inversion, these bounds mean in real time 

C 3 t(logt) 1/2 < E (| X(t) | 2 ) < C 4 t logt, 
with < C3, C4 < oo and for t sufficiently large. 

2 Fock space computations 
2.1 Spaces and operators 

For basics of Gaussian Hilbert spaces (that is: Fock spaces) see [17] or |7J. Let (Q, 7r) be 
either one of the probability spaces identified in subsection 11.31 The Gaussian Hilbert 
space "H = C 2 (Q, tt) is naturally graded 

oo 

H — Tin 

n=0 

where 

Ti n = span{:w ;i (xi) . ..ui n (x n ): : h, ...,/„ = 1, 2, xx, ...,x n e R 2 }. 

Here and in the sequel : Z\ . . . Z n : denotes the Wick monomial formed by the jointly 
Gaussian random variables Z±, . . . , Z n . 

The operators k = 1,2, and A are the usual ones. Their action on Wick mono- 
mials is 

n 

V fc :uj h (xi) . . . u ln (x n ): = ^ :UJ h (Vi) ■ ■ ■ d ^i m (x m ) ...u ln (y n ): , 

m=l 

A = IX 

k=l 

They are defined as unbounded operators on Ti by graph closure. We will also need 
the following creation and annihilation operators. Again, we give their action on Wick 
monomials: 

a* k :u h (xi) . ..ui n (x n ): =:uj k (ti)u h (xi) . ..ui n (x n ): 

n 

a k :u h (x l )...UJi n (x n ):= ^K klm {x m ) w h (y x ) . . . u^J^x^j . . . uj ln (y n ): . 

m=l 

Note that the creation and annihilation operators defined in the context of the different 
models models are not unitary equivalent. 
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In both isotropic cases (SRBP and DCGF) the infinitesimal generators of the semi- 
groups of the stationary Markov process t] t , acting on £ 2 (Q, ir), is expressed as 

G = -S + A + + A_, 

2 2 

S = -A, A + = 4V fc , A - = Y1 Vfcflfc - 

k=l k=l 

This follows from standard computations in the case of DCGF. For the SRBP process 
it relies on somewhat more complex considerations with some involvement of Malliavin 
calculus. For details see [6]. Note, that the infinitesimal generators of the two processes, 
although formally similarly expressed, are not unitary equivalent. 

For the non-isotropic SRBP, since w^ipc) = 0, a\ = = 0, the asymmetric part of 
the generator is slightly modified: 

A + = a*Vi, A_ = Viai. 



2.2 The variational formula 

We write the displacement as sum of a martingale and compensator term. In each of the 
cases we will concentrate on the first coordinate component. 

X l (t) = B l (t)+ [ <p( Vs )ds, (7) 

where 

Note, that actually tp G Tti. 

The martingale term is diffusive, so in order to prove superdiffusive bounds we have 
to focus on the integral on the right hand side of fl7]). Laplace transformation yields: 



rt \ 2 N 

A< E ( ; / tp{r)s)ds) ) dt = \~ 2 (<p 
where 

R x = e~ xt e tG = (XI-G)-\ 
Jo 

is the resolvent of the infinitesimal generator G. We are going to prove bounds for the 
right hand side of (jSJ). 

The following variational formula is straightforward, see e.g. [12J: 

(<p, R xV ) = sup {2(<p, 4>) - (V, (A + 5)^) - (Aip, (A + S)- l A^) } 

The upper bounds in (jSJ), © drop out essentially for free: 

(<p, R xV ) < sup {2(<p, ^) - (i>, (A + S»} = (<p, (A/ - A)"V) < Clog A. 



The last bound follows from straightforward computations on so-called diffusion in ran- 
dom scenery, and proves the upper bound in both Theorem [T] and [21 
In order to get the more interesting lower bound we write first 

{ip,R*p)> sup {2(^^)-(^(A + 5)^)-W,(A + 5)- 1 ^)} 
i/>eHi 

= sup {2( l p^)-(^(X + S)iP)-(A + iP,(X + S)- 1 A + iP)}. (9) 

Note that for any ip G "Hi, A^tp = 0. 

We will treat the isotropic SRBP and DCGF processes and the anisotropic SRBP 
separately. 

Isotropic SRBP and DCGF models: 

We write the variational test function ip G Hi as 

2 „ 

ui(x)ui(x)dx 



< oo, 



*(<»>) = E / 

where u : IR 2 — > M? is such that, u(—x) = u(x), (and thus u : M 2 — > M 2 ), and 

SRBP : p x u(p) = 0, / V(p) (p ■ u{p)f dp 

Jr. 2 

DCGF: p-u(p) = 0, / V(p) (p x u(p)fdp < oo. 

We define u : M 2 -> i R be 

SRBP: t> = idivw, v{p)=p-u{p), 

DCGF: -y = zrotu, u(p) = p x u(p). 

Then, in both cases v : R 2 — )■ IR and 

v(—p) = —v(p), and / jv(p) 2 dp < oo. (10) 

Jr 2 \p\ 

The first two terms on the right hand side of the variational lower bound Q are 
computed directly. 

SRBP: M)= ( V(p)-^v(p)dp=:J 1 (v) (11) 

Jr 2 \p\ 

DCGF: M)= ! V(p)-^v(p)dp=:J 1 (v) (12) 

Jr. 2 \p\ 



SRBP&DCGF: (i/>,(\ + S)i/>)= V{p) ^ v(p) 2 dp =: J 2 (v). (13) 

Jr 2 \p\ 
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In order to compute the third term on the right hand side of §§§ we first consider the 
SRBP case. We have 



- y~] / d k ui(x) :u k (0)ui{x): dx, 
k,i=i jR2 

and hence 

(A + ij,(\ + Sy 1 A + ^)= (14) 
/ / / d m u n (y)d k ui(x)g\(z)E(:u m (0)u n (y)\\u k (z)ui(x + z)\) dxdydz. 

Here g\(z) is the integral kernel of the operator (A — A) -1 in IR 2 , with Fourier transform 

9*(p) = x , \ |2 - 
A+ \p\ 

For the DCGF case one finds very similar formulas, with d k replaced by d k . 

The expectation inside the integral in (|14|) is computed by exploiting the fact that the 
fields u are Gaussian and hence the four point functions arising are expressed in terms 
of the covariances. After some straightforward computations, eventually we get: 

SRBP: {\ + S)- 1 A + il)) = (15) 

V{p)V{q) (P 2 f 1 - 2 (v(p) - v{q)fdqdp =: J 3 (v), 

\p\ \q\ A + \ p — q\ 

DCGF : (A+i/>, (A + S^A+^p) = 

V(p)V(q) ( P * q)2 1 ,(0(p) -%)) 2 ^P=: J 8 («). 

\p\ \q\ a + \ p — q\ 

And the variational problem becomes 

(<p,Rx<p) > sup (2J x (v) - J 2 (t)) - J 3 (v)), (16) 

i 

with supremum taken over functions v : M 2 — )■ R satisfying conditions ( ITU]) . 

Given the explicit expressions (fTTj) . ([121) and ([13]) for Ji('O) and J2{v), in order to 
prove the superdiffusive lower bound in we need to prove efficient upper bounds on 

Hv). 

Anisotropic SBRP case: 

The test function i/j e Tii is written in the form 



i/j(u:) = / u(x)ut(x)dx 
11 



with u : R 2 -> R, m(x) = u(-x). Then 

{f^)= I u(p)V(p)dp =: Ji(w), 
Jr 2 

(^,(A + £» = / V(p)(A+|p| 2 )fi(p) 2 dp=: J 2 (m) 
Jr 2 

(Ah^(A + 5) _1 ^)= / / \>(p)\>(g)( Pl n(p)-g 1 n(g)) 2 ^g=: J 3 (u). 

7r 2 Jr 2 

3 Computations 

The following three subsections contain the computational parts of the proofs for the 
isotropic cases (subsections 13.11 and 13. 2p . respectively the anisotropic case (subsection 

ED- 

3.1 Upper bound on J3, isotropic models 

Since the DCGF model is somewhat simpler than the SRBP, we will treat them in this 
order. 



DCGF: 

Denote 



D(A,|p|)=4 / V(q) 



(p x q) 



2 1 

2 1 72 7 ~\ |2 ( 



I P I I 9 I A + |p — g 

When estimating the function D(X, \p |) (see Lemma [T] below and its proof) the vectorial 
product p x q will help substantially. As we shall see later this is not the case for the 
SRBP model. 

By Schwarz's inequality we get 

Mv) < [ V(p)D(\,\p\)v(p) 2 dp. 

and the variational problem (TTB]) is readily solved by 

v \P) — o, 

A+(1 + D(A,|j>|))IpI 
which eventually yields the lower bound 

(<p,Rx<p)> [ V(p)(X + (l + D(X,\p\))\p\ 2 y 1 dp. (17) 

Thus, a suitable upper bound on D(X,\p\) will provide a good lower bound on (ip, R\<p). 
In subsection 13.21 we will prove the following 
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Lemma 1. In the DCGF case, for A < 1 and \ p \ < 1 we have 

D(X,\p\)<C\\og(X + \p\ 2 )\. (18) 
Inserting (lT8"j) into (TT7|) we readily obtain 

(<p, R x <p) > [ l iM<1} V(p) (A + (1 + C | log(A + | p | 2 ) | ) | p | 2 ) -'dp > C log | log A | . 
Jr 2 

which proves the lower bound in Theorem [1] for the DCGF. 
SRBP: 

Applying a similar argument to the SRBP process will not yield super diffusive lower 
bound: the bound corresponding to ffT8|) would be 

D(A,b|)<C|logA| 

and this is simply not sufficient for superdiffusivity. 

In order to get the true super diffusive lower bound we have to split the integral on the 
right hand side of f)15p according whether | p — q | is small or large and apply Schwarz's 
inequality only in the latter part. Let 

Jsi(v) ■■= I I V(p)V(q) #(v(p) - v(q)) 2 l {lp ^ ql > lpl/3} dqdp 

7r2 7r2 A + | p — q | 

Js2(v)= I I V(p)V(q) (v^p) - v(q)) 2 l {lp _ q \< lpl/3} dqdp 

Jw? Jm? A + | p — q | 

Then, by applying Schwarz's inequality 

Mv)<2 [ [ V(p)V(q)- p ^{v{pf + v{qf)t { \ p ^< Vpm dqdp 

< f V(p)D(\,\p\)v(p) 2 d Pj 
Jr 2 

where now 

D(X,\p\) = A / V(q)—— ^{\p- q \>\ P \/3}dq. 

Jr 2 A + | p — q I 

In subsection 13.21 we prove the following 

Lemma 2. In the isotropic SRBP case, for A < 1 and \ p | < 1 we have 

D(X,\p\) <C|log(A+M 2 /9)|. 
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Now, we give upper bound on J 32 : 

Js2< I I V(p)V(q)- — -2 ((p- q) ■ Vv(r(p,q))) 2 l { \ p - q \<\ p \ /3} dqdp 

7r2 J R 2 A + I p — q | 

- 7 / ' P I 2 SUP I V ^( r ) I 2 dp =: J 32(^) 

4 </IR 2 r:\r-p \<\p\/3 

We choose the variational function 

v*(p) := cpih(X + |p| 2 ), h(x) := ,\ (19) 

Then at one hand, for c sufficiently small we get 

Ji(v*) ~ ~ Jsi{v*) > C log I log A I . (20) 

On the other hand, in the next subsection we also prove 
Lemma 3. With the choice ()19j) of the variational function we have 

4 2 (v*) < C. (21) 
Finally, from (J2TJ and ([21]) it follows that 

(<p, R x <p) > C log | log A | 
which gives the lower bound Theorem [T] for the SRBP. 

3.2 Proof of Lemma Q} El and G2 

Proof of Lemma [7J By the rotational symmetry we may assume that p > is real. Using 
the decay of V the integral on \q\ > 2 is bounded by a fixed constant. Since V is bounded 
it is enough to bound 

f (p x q) 2 i 



2 , , 72^- (22) 

<2 | p | | 9 | A + | p - g I 



Let q = (r cos(t), r sin(i)) then we may rewrite (|22[) as 

,2 .2y 2 in 2M j ,2 ,2. Sin 2 (t) , , , . 

/ / ?rdtdr = / / - —rdtdr. 23 

Jo Jo P r X+\p-re u \ Jo Jo A + r 2 + p 2 - 2pr cos(t) 

If A > B > then 

sin 2 (t) sin 2 (7r + t) 2 A sin 2 (i) 2 

A + B cos(t) + A + 5 cos(vr + t) ~ A 2 - 5 2 cos 2 (t) ~ A 

which shows that after integrating in t in (l2"3"j) : 

r 2 r 2 " sm 2 m , , r 2 2tt , /A 2 ^ , /y 

Collecting all our estimates and using A < 1, |p| < 1 the statement of the lemma follows. 

□ 
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Proof of Lemma\^ By the decay of V the integral on \p — q\ > 2 is bounded by a fixed 
constant. Thus we only need to bound 

1 f 2 1 

V ^V~\ ^{\P-i\>\ P m d( l = / YT^ xdx - 

<\p-q\<2 X+\p — q\ J\p\/3 A + X 

The last integral is log (A + 2) — log (A + |p| 2 /9) from which the lemma follows. □ 

Proof of Lemma^ Recall the definition of the variational function v* from (TIT)]) . In order 
to avoid heavy notation we will drop the * from v* in the subsequent computations. We 
have 

Vv(p) = (h(X + \p\ 2 ) + 2p\h\X + | p | 2 ), 2p 1 p 2 h\X + | p | 2 )) 

and hence 



I Vv(p) | 2 < 2h(X + | p | 2 ) 2 + 8\p\ 4 h'(X + | p | 2 ) 2 

< 2h(X+ |p| 2 ) 2 + 8(A + \p\ 2 ) 2 h'(X + \p\ 2 f. 

From the explicit form of h: 

1-x 2 



(24) 



h'(x)x 



1 



< 2 



x(l + x 2 ) log 2 (- + x) x log + x) 
which together with (12~4"]) yields 

I Vt)(p) | 2 < 34/i(A+ |p| 2 ) 2 

Thus it is enough to bound 



1 



xlog (| + x) 



= 2\h(x) 



V(p) \p\ 2 sup 

r:| r—p \<\p |/3 



/i(A + | r | 2 ) 2 (ip. 



Note that from the definition of h(x) (see (I19p ) we have 

,2x2 



| p | sup h(X + I ■ 

r:| r—p |<| p |/3 



< 81- 



p| 2 log 2 (|p| 2 /9 + 9| P r 2 )' 



(25) 



The lemma now follows from the fact that V(p) is bounded and the function on the right 
hand side (125 j) is locally integrable near | p | = in M 2 . □ 



3.3 The non-isotropic SRBP 

The proof of the lower bound is considerably simpler in that case. We use the same 
strategy as before, setting 



D(A,|p|) = 4 / V(p) 



1 



A+ \p — q\ 2 



dq 
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by Schwarz's inequality we get 

Js(u) < f V(p)D(\,\p\) P lu(p) 2 dp. 
Jr 2 

As in the DCGF case we have that if A < 1 and \p\ < 1 then 

D(X, M)<C|log(A+p 2 )|. 

Using the same arguments as before we get the following lower bound for sufficiently 
small A: 



(<P,R\<p) > / *{\ P \<i}V(p) — — — ——r 2 . dp. 

A + \p\ 2 + Cp{ log (A + | p | ) 



1 

A + |p| 2 + Cpf |log(A) 



> / l{|p|<l/2}V"(y) - ; ; 7^2TTZZ7m d P- 



Changing to polar coordinates and using that V(p) > C > for \p\ small enough we 
have the lower bound 

e r2n r , . r 2ixr 



C' I I dadr = C' dr 

J Jo Jo A + r 2 + C|logA|r 2 cos 2 « J J ^(r 2 + A) (C\ log A|r 2 + r 2 + A) 

> C'^loiA]. 

This proves Theorem [21 

4 Appendix: The Alder- Wainwright scaling argument 
for superdiffusivity of tracer motion 

We reproduce the nonrigorous, nevertheless very instructive scaling argument due to B.J. 
Alder and T.E. Wainwright, respectively, to D. Forster, D. Nelson and M. Stephen, which 
sheds sharp light on the origins of superdiffusivity in tracer particle motions. The original 
papers are PQ, j2], [5], see also [18J. 

4.1 General formal setup and notation 

Let 1 1— > X(t) G K d be a random motion of a tracer particle, with stationary and ergodic 
increments. The motion is performed in some random environment which also evolves in 
time. The time evolution of the tracer particle and that of the background environment 
may mutually influence one another. 

Usually we decompose the random motion as sum of a martingale and its compensator: 

X(t) = M{t) + f V(s)ds, (26) 
Jo 
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where M(t) is a square integrable martingale with stationary and ergodic increments and 
the compensator V(t) is a stationary and ergodic process which is also square integrable. 
We'll call V(t) the instantaneous velocity of the tracer. We are interested in understand- 
ing the origins of the possibly superdiffusive behaviour of X(t). Since M(t) is anyway 
diffusive, we don't care much about it. The main contribution to the superdiffusive be- 
haviour anyway comes from the compensator (integral term) on the right hand side of 



It is assumed that the instantaneous velocity comes from some background velocity 
field U (t, x) as follows: 

V(t) = U(t,X(t)). (27) 

It is important to note that we think about a joint random dynamics t h- >• (X(t), U (t, •)) . 
The process of the velocity field t (->• U(t, •) is usually not stationary while the instanta- 
neous velocity of the tracer given in (|27|) is. 

We will consider here only the isotropic case. The results of similar considerations 
applied to non-isotropic cases will be also summarized at the end of this appendix. 

The correlations of the velocity field: 

K(t,x) : = E (£7(0,0) • U(t,x)) . 
The velocity autocorrelation function: 

C{t) :=E(V(0)-V(t)) 
Mind that the velocity process V(t) is assumed stationary (and ergodic), thus 

E(V(s)-V(t)) =C(t-s). 
The variance of the displacement is: 



E(t) := E\X(t) \ 2 , E(t):=E 



V(s) 



Since the martingale part in (1261) is anyway diffusive, in case of superdiffusive behaviour 

E(t) x E(t). 

By stationarity of V(t) 

ft 



E(t) = 2 I (t- s)C{s)ds, 
Jo 



and hence 



E(t) x t I C(s)ds. (21 



We shall refer to (I25|) as the asymptotic form of the Green-Kubo formula. 
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4.2 Scaling assumptions and computation 
Scaling of the displacement: We assume that X(t) is of order 

a{t) =r(logt) 7 . 

More explicitly, assume that 

P (X(t) E dx) x a(t)- d ip{a(t)' 1 x)dx, (29) 

where </? : M. d — > M + is a density, which is regular at x = and decays fast at | x \ — > oo. 

Scaling of the velocity field: We also assume that the correlations of the velocity field 
U(t,x) scale as 

K(t,x)-m- d Hm~^)- 

Note that under this assumption 

/ K(t, x)dx x const. 

This corresponds to some kind of conservation of momentum carried by the velocity field. 
Main scaling assumptions: We make two important assumptions 

1. Regularity of and tp: The following regularity conditions hold: 

\^(x)\dx<oo, V>(0) > 0. (30) 



2. Displacement scales on faster order than the velocity-field correlations: 

/?(*) = 0(a(t)). (31) 

Next we compute the velocity autocorrelation function: 
C(t) = E 0/(0) • V(t)) = E (17(0, 0) • U(t, X(t))) 

/ K(t,x)~P (X(t) (E dx) decoupling 

/3(t)- d ip(/3(t)- 1 x)a(t)- d (p(a(t)- 1 x)dx scaling 

= a{ty d I ^j(x)y(f3{t)a{ty 1 x)dx 

x a(t)~ d (32) 

Of course, the second step (decoupling) is the shaky one. In the last step we used the 
main scaling asumptions ( 13T)]) and (13"T|) . Regularity of cp at x = was also assumed. 

18 



/ 



4.3 Conclusions 

From the scaling assumption f l29|) it follows that 

E(t) x a(t) 2 . (33) 



On the other hand, using the Green-Kubo formula (128]) and the computations in fl32|) we 
get 



t 

d. 



E(t) x t j a(s)- d ds. (34) 
The only choices of the scaling function a(t) consistent with both fl33l) and (|34|) . are 

d=l: u=l, 7 = 0, £(t)xt 4 / 3 , 

d = 2: 7 = i S(t)x*(log0 1/2 , 

d = 3: !/ = - 7 = 0, E(t)xt. 

If d = 2 and the system is non-isotropic in the sense that X^it) is assumed diffusive 
and Xi(t) possibly superdiffusive, very similar considerations and computations lead to 
the following conclusion: 

d = 2: „=1 7 = I E{t)^t{\ogt) 2 '\ 
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